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Abstract 

We study the quantum paraelectric-ferroelectric transition near a quantum critical point, em- 
phasizing the role of temperature as a "finite size effect" in time. The influence of temperature 
near quantum criticality may thus be likened to a temporal Casimir effect. The resulting finite- 
size scaling approach yields ^ behavior of the paraelectric susceptibility (x) and the scaling form 
x(w,T) = -^jF{^), recovering results previously found by more technical methods. We use a 
Gaussian theory to illustrate how these temperature-dependences emerge from a microscopic ap- 
proach; we characterize the classical-quantum crossover in x, and the resulting phase diagram is 
presented. We also show that coupling to an acoustic phonon at low temperatures (T) is relevant 
and influences the transition line, possibly resulting in a reentrant quantum ferroelectric phase. 
Observable consequences of our approach for measurements on specific paraelectric materials at 
low temperatures are discussed. 
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I. INTRODUCTION 



The role of temperature in the vicinity of a quantum phase transition is distinct from 
that close to its classical counterpart, where it acts as a tuning parameter. Near a quantum 
critical point (QCP), temperature provides a low energy cut-off for quantum fluctuations; 
the associated finite time-scale is defined through the uncertainty relation At ~ j^t- This 
same phenomenon manifests itself as a boundary condition in the Feynman path integral; 
it is in this sense that temperature plays the role of a Unite-size effect in time at a quantum 
critical point .-^^ The interplay between the scale-invariant quantum critical fluctuations 
and the temporal boundary condition imposed by temperature is reminiscent of the Casimir 
effect,-^ where neutral metallic structures attract each other-^^^^ due to zero-point 
vacuum fluctuations. 

In this paper we explore the observable ramifications of temperature as a temporal 
Casimir effect, applying it to the example of a quantum ferroelectric critical point (QFCP) 
where detailed interplay between theory and experiment is possible below, at and above the 
upper critical dimension. Our work is motivated by recent experiments on the quantum para- 
electric SrTiO-j (STO) where 1/T 2 behavior is measured in the dielectric susceptibility near 
the QFCP.— dkdS. Here we show how this result is simply obtained using finite-size scaling in 
time; more generally we present similar derivations of several measurable quantities, recov- 
ering results that have been previously derived using more technical diagrammatic,— 
large N— and renormalization group methods.— 1 ^ In particular we present a simple inter- 
pretation of finite-temperature crossover functions near quantum critical points previously 
found using e-expansion techniques,— and link them to ongoing low-temperature experi- 
ments on quantum paraelectric materials. We illustrate these ideas using a Gaussian theory 
to characterize the domain of influence of the QFCP and we present the full phase diagram. 
Next we expand upon previous work by tuning away from the QFCP, studying deviations 
from scaling; here we find that coupling between the soft polarization and long- wavelength 
acoustic phonon modes is relevant and can lead to a shift of phase boundaries and to a 
reentrant quantum ferroelectric phase. We end with a discussion of our results and with 
questions to be pursued in future work. 
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II. THE CASIMIR EFFECT IN SPACE AND IN TIME 



The Casimir effect is a boundary condition response of the electromagnetic vacuum. The 
gapless nature of the photon spectrum means that the zero-point electromagnetic fluctu- 
ations are scale- invariant; the vacuum is literally in a quantum critical state. However, 
once the boundary conditions are introduced, the system is tuned away from criticality and 
develops a finite correlation length, £. The Coulomb interaction between two charges, the 
correlation function of the electromagnetic potential inside the cavity, is changed from the 
vacuum to the cavity as 

V(q)f r ee ~ (50 q 50- q ) = — -> V{q)cavity 2 ^2 f = ~ ( 1 ) 

where the plates have removed field modes and have introduced a finite £. In an analogous 
way, the partition function of a quantum system at finite temperatures is described by a 
Feynman path integral over the configurations of the fields in Euclidean space-time^ where 
temperature introduces a cutoff in the temporal direction. In Figure 1 we present a visual 
comparison of the Casimir effect in space and in time. In both cases, the finite boundary 
effects induce the replacement of a continuum of quantum mechanical modes by a discrete 
spectrum of excitations. 

In the quantum paraelectric of interest here, the path integral is taken over the space-time 
configurations of the polarization field P(x, r), 

S E [P] 



5^ exp 

{P(x,r)} 



h 



(2) 



where 

S E [P] = / kBT drd 3 xC E [P} (3) 
Jo 

and £e[P] is the Lagrangian in Euclidean space-time. The action per unit time is now the 
Free energy F of the system (See Table I.). The salient point is that finite temperature im- 
poses a boundary condition in imaginary time and the allowed configurations of the bosonic 
quantum fields are periodic in the imaginary time interval r G [0, h/3] ((3 = ■j^f) so that 
P(x, t) = P(x, t + h/3), which permits the quantum fields are thus decomposed in terms of 
a discrete set of Fourier modes 

P n (x, t)=J2 ^n)e iM (4) 
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Figure 1: Casimir effect in space and time, (a) Imposition of spatial boundaries on the quantum 
critical electromagnetic field yields the conventional Casimir effect (b) Imposition of temporal 
boundary conditions on a quantum critical paraelectric generates the effect of non-zero temperature. 



where 



v n = n 



2irk B T 



(5) 



are the discrete Matsubara frequencies; we recall that at T = the (imaginary) frequency 
spectrum is a continuum. The response and correlation functions in (discrete) imaginary 
frequency 



XE(q,iv n ) = (P{q,iv n )P{-q,-iVn)) 
can be analytically continued to yield the retarded response function 



(6) 



XE{q,iVn) -> XE(q,u) = Xe(<1, Z)\z=u+i6 



(7) 



where u is a real frequency; for writing convenience we will subsequently drop the "E" 
subscript in x e.g. x(<7» = 

Table. 1. Casimir Effect and Quantum Criticality. 
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Like the parallel plates in the traditional Casimir effect, temperature removes modes 
of the field. In this case it is the frequencies not the wavevectors that assume a discrete 
character, namely 

q = (q,u) -»• (<f,iv n ), (8) 



where v n are defined in (JSJ). 

The Casimir analogy must be used with care. In contrast to the noninteracting nature 
of the low-energy electromagnetic field, the modes at a typical QCP are interacting. In 
the conventional Casimir effect, the finite correlation length is induced purely through the 
discretization of momenta perpendicular to the plates. By contrast, at an interacting QCP, 
the discretization of Matsubara frequencies imposed by the boundary condition generates 
the thermal fluctuations in the fields in real time. These are fed back via interactions 
to generate a temperature-dependent gap in the spectrum and a finite correlation time. 
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Despite the complicated nature of this feedback, provided the underlying system is critical, 
temperature acting as a boundary condition in time will set the scale of the finite correlation 
time 

tr = -T^, (9) 

where k is a constant. In cases where the quantum critical physics is universal, such as 
ferroelectrics in dimensions below d = 3, we expect the coefficient k to be also universal 
and independent of the underlying strength of the mode-mode coupling. The "temporal 
confinement" of the fields in imaginary time thus manifests itself as a finite response time 
in the real-time correlation and response functions. 

For the quantum paraelectric at the QFCP, the imaginary time correlation functions are 
scale-invariant 

x (q,iv)= (P(q,iv)P(-q,-iv)) 



~ 9 1 9 9 - (10) 



At a finite temperature this response function acquires a finite correlation time 

1 

v 2 n + c 2 s q 2 + Cr 



where 

C 2 = 31c{(P 2 )t^-(P 2 )t=o} (12) 

is determined by mode-mode interactions, where 7 C is the coupling constant describing the 
quartic interactions between the modes, to be defined in Sec IV. We note, as shall be shown 
explicitly in Section IV, that for dimensions d such that 1 < d < 3, the feedback will be 
sufficiently strong such that £ T will be independent of the coupling constant 7 C ; by contrast 
for d > 3 the feedback effects are weak so that there will be a 7 c -dependence of £ T . The case 
d — d% — 3 is marginal and will be discussed as a distinct case. At a temperature above a 
quantum critical point, the energy scale 

A(T) = ak B T (13) 

will set the size of the gap in the phonon dispersion relation. Here A(T) ~ 1 and 
a = 0(1) is a constant of proportionality. 

Reconnecting to our previous discussion, we remark that real-time response functions 
from expressions like (fTTj) are obtained by analytic continuation to real frequencies %v n — > uj. 
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q=0, v =0 



in contrast to the Curie form (x ~ t) associated with a classical paraelectric; this 1/T 2 
temperature-dependence was previously derived from a diagrammatic resummation,— 
from analysis of the quantum spherical model 2 ^ and from renormalization-group studies.— 1 ^ 
We note that this 1/T 2 behavior in the dielectric susceptibility of the quantum paraelectric 
has been observed experimentally.— 1 ^ 1 ^ We summarize in Table I the link between the 
conventional Casimir effect and finite-temperature behavior in the vicinity of a QCP. 

III. FINITE-SIZE SCALING IN TIME 

The spatial confinment of order parameter fluctuations near a classical critical point 
has been studied as a "statistical mechanical Casimir effect" and here we extend this 

treatment to study the influence of temperature near a QCP using finite-size scaling (FSS) 
in imaginary time. This scaling approach is strictly valid in dimensions less than the upper 
critical dimension. Quantum critical ferroelectrics in d = 3 lie at the marginal dimension 
(D = d + z = 4), so the scaling results are valid up to logarithmic corrections, which we 
discuss later (Sec. VI); here z — 1 refers to a linear dispersion relation, to = c s q. 

Following the standard FSS procedure,- 1 ^ 1 ^ we impose boundaries on the system near 
its critical point. For a classical system with tuning parameter t = T ~ Tc and correlation 
length £, we confine it in a box of size L and then write the standard FSS scaling form 



Here temperature is no longer a tuning parameter, this role is taken over by an external 
tuning field g. Temperature now assumes a new role as a boundary condition in time. 
Introducing a fixed L T (see Fig. 2b) associated with a finite T, while replacing t —* g, the 
quantum critical version of (115j) is 




(15) 



for the susceptibility 



Similar reasoning can be used when a system is near its QCP. 




(16) 
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b) 




L T » £ T 



Figure 2: Schematic of finite-size effects (a) at a classical and at a (b) quantum critical point where 
the appropriate lengths are defined in the text. 

where g is the tuning parameter. The dispersion relation u> = c s q z yields [£ T ] = [£*]; this 
combined with £ ~ g~ v leads to £ T ~ g~ zu ■ We therefore write 



where <$>(x) ~ x p is a crossover function where p is determined by the limiting values of 
when x — > 0, we expect x = x(-^r); whereas we should recover the zero-temperature 
result (x ~ S 1 " 7 ) when x — > oo. Therefore we obtain 



and the temperature-dependence (L T oc 1 /T) emerges naturally from FSS arguments. There- 
fore a (T = 0) quantum critical point can influence thermodynamic properties of a system 
at finite T just as a finite-size system displays aspects of classical critical phenomena despite 
its spatial constraints. A schematic overview of the finite-size scaling arguments we have 
presented here is displayed in Figure 3. 

The FSS approach can also yield the T-dependences of the specific heat and the polariza- 
tion of a quantum critical paraelectric. At a finite temperature phase transition, to obtain 




(17) 




(18) 
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Figure 3: Overview of the finite-size scaling at classical and quantum critical points; here i> is the 
exponent associated with the spatial correlation length since v has already been used in the text 
as a frequency. 

the specific heat capacity of a finite size box with L « £, we write / ~ t 2 ~ a F(L/C,) ~ 
t 2 ~ a (prj) 2 ~ L-( 2 - a )/ u . In a similar spirit, applying the quantum FSS analogies 
(L -> L r , i -»• £ ^ = #~ 2p ), we obtain 

fULr)-9 2 - a {^j " -L^-T^ (19) 

so that the T-dependent specific heat is 

c v (T)=T&f qm /8Z a ~T s &- 1 (20) 

in the vicinity of a QCP. Similarly the T-dependence of the polarization is P{T) ~ 
and we note that P(E) = df qm /dE\ g=0 ~ E 1 / 5 is T-independent, since finite-temperature 
scaling does not affect field-behavior. 

Simple scaling relations at classical and quantum criticality are summarized in Figure 3. 
The key notion is that at a QCP, finite T effects correspond to the limit L T « £ r ; in this 
case L T becomes the effective correlation length in time, and the T-dependences follow. We 

9 



note that we expect this finite-size approach to work for dimensions d < d u c where there will 
be logarithmic corrections to scaling in the upper critical dimension d u c . 

Let us now be more specific with exponents for the quantum paraelectric case. At crit- 
icality the observed T-dependence of the paraelectric susceptibility (x) can be found by a 
soft-mode analysis,— ^ and therefore the exponents for the quantum paraelectric are those 
of the quantum spherical model.— For the case of interest (D = d + z = 3 + l = 4), the 
quantum spherical model has exponents v = 1/(D — 2) = 1/2 and 7 = 2/(D — 2) = 1, 
so that = 4r = 2 and we recover the \~ x ~ T 2 scaling found earlier. Other specific 
T-dependences are displayed in Table II. For d — 3, we have g ~ T 2 ; this relation was 
experimentaly observed^ 3 - in STO. Finally we note that the FSS that we have discussed 
suggests the "—" scaling form 

that is similar to that observed in other systems at quantum criticality;— ^ this was previ- 
ously derived by more technical methods.— Predictions for experiment are summarized in 
Table II. We note that since we are in the upper critical dimension, there will be logarithmic 
corrections to this scaling but we do not expect these to be experimentally important for 
the temperature dependences described here; however they will be considered later in the 
paper (Section VI). 

Table II. Observables for a D = 4 QPE in the Vicinity of a QFCP 



Observable 


T-Dependences 

(g=o) 


g-Dependences 
(T=0) 


Polarization 


P-T 1 


P ~ gh 


Susceptibility 


X~T- 2 


x ~ g' 1 



p ~ eI 

X (^,T) = ^F(|) 
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IV. GAUSSIAN THEORY: ILLUSTRATION OF TEMPERATURE AS A BOUND- 
ARY EFFECT 



A. The Gap Equation 

In this section we use the self-consistent Gaussian theory to illustrate how the x(T) found 
via FSS in time appears from a more microscopic approach; we also study the crossover be- 
havior between the classical and the quantum critical points. This approach is equivalent 
to the self-consistent one-loop approximation^ that is used in the context of metallic mag- 
netism. 

The soft-mode treatment has been described extensively elsewhere;— i22a2& here we briefly 
outline the derivation of the gap equation. The Lagrangian in Euclidean space-time, Ce in 
dSJ), for displacive ferroelectrics is the 4 model: 

- \ [(d T P) 2 + (VP) 2 + rP 2 ] + |P 4 . (22) 

which determines the partition function. Notice that in writing (1221) . we have chosen rescaled 
units in which the characteristic speed of the soft mode c s — 1. In a self-consistent Hartree 
theory, interaction feedback is introduced via its renormalization of quadratic terms; this 
procedure is equivalent to replacing Ce by the Gaussian Lagrangian 

C G = l -P[-d 2 T -V 2 + r + Y]P (23) 

where 

£ = 3 7c (P 2 ) (24) 

is the Hartree self-energy (see Fig. 4). We note that this mode- mode coupling theory is 
exact for the "spherical model" generalization of 4 theory in which the order parameter 
has N components and N is taken to infinity. 

The Green's function can now be determined from Dyson's equation, shown diagrammat- 
ically in Figure 4, and takes the form 

G(q) = G(q, iv n ) = [(iv n ) 2 - g 2 - r - S] - 1 , (25) 

so the action is diagonalized in this basis. The poles of G(q, uS) determine the dispersion 
relation uj q for the displacive polarization modes 

uj\ = q 2 + A 2 (26) 
11 
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Figure 4: Diagrammatic Representation of (a) the Dyson Equation and (b) the Gaussian Self- 
Energy where Pq = in the paraelectric state (and is finite in the ferrelectric phase. 

where here we have introduced the gap function 

A 2 (r,T) = r + £(r,T). (27) 

This quantity vanishes at both quantum and classical critical points where there are scale- free 
(gapless) fluctuations. At the quantum critical point where T c = 0, A(r c , 0) = r c + £(r c , T = 
0), so that we can eliminate r c = — £(r c ,T = 0), to obtain 

A 2 (r, T) = Q 2 + [E(r, T) - £(r c , 0)], (28) 

where Qq = (r — r c ) = g. 

The amplitude of the polarization fluctuations is given by 

<P 2 > = -G(0, 0) = --L ^ G(g)e^U =0 , (29) 
so the self-consistency (1241) condition S = 3 , y c {P 2 ) can now be written 

S(r,T) = (" 3 7c)^E / (^ G ^ ZZ/ «)' ( 3 °) 

where S(r, T) is the temperature-dependent self-energy. By converting the discrete Matsub- 
ara summation to a contour integral, deformed around the poles z = ±ui(q) in the dispersion 
relation, we can convert this expression to the form 
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where we denote ub{oj) = ^(w,/?) = l/(e l3uJ — 1). At the quantum critical point (r = r c 
and T = 0), we have u q = q and n(uj q ) = so that 

and using (1281) . we can write the gap function as 



rc 7 (2tt)^ I 2g 



V 7 ? 2 + A 2 . (33) 



B. T-Dependence of the Gap at the QCP. 



In the paraelectric phases, we can use the temperature-dependent gap to determine the 
dielectric susceptibility x- Writing 



[P(q)P(-q)) 



q,u)=0 



q=0 



(34) 



q,cv=0 



we use ( 1251) and (1271) to express it as 



(35) 



X ~\r,T) = A 2 (r,T). 
At the quantum critical point SIq = 0, so the gap equation (I3"3"j) becomes 

Jo (2ir)*l ^ + 2ql' 1 ' 

where we have inputted the dispersion relation, (126]) . for uo q in (136]) . We notice that both 
thermal and quantum fluctuations contribute to this expression. 

Even though the mean field gap equation is only formally exact in the spherical mean- 
field limit, it is sufficient to illustrate the qualitative influence of T on the gap at the QCP. 
In order to explore the cutoff-dependence of (1361) . we note that in the ultraviolet limit of 
interest, the last two terms can be expressed as 

\{k~ l ~q) = ~W (37) 

where there is complete cancellation when A = exactly at the QCP. However just slightly 
away from it, when A is finite, (l3Tj) leads to a qf^x scaling-dependence of the integral in (1361) ; 
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therefore the cutoff is required to ensure that (136]) is finite in dimensions d > 3. However, in 
dimensions d < 3, this integral is convergent in the ultraviolet and the upper cutoff in (136]) 
can be entirely removed. Thus, for d < 3, the only scale in the problem is temperature itself. 
The integral is also convergent in the infrared provided d > 1 . The spatial dimensions d — 1 
and d = 3 correspond to spacetime dimensions D = d + 1 = 2 and D = d + 1 = 4, which 
are the well-known lower and upper critical dimensions of the 4 theory. This provides us 
with a dimensional window 1 < d < 3 where inverse temperature acts as a cut-off in time. 
In this range, the temperature-dependence of the gap 

A(T) = a d T (38) 

is independent of the strength of the coupling constant 7 C and the cutoff, a feature that 
can be illustrated already within mode-coupling theory. Recalling that A(T) = -f- and 
a = aa (see ([13]) and Fig.l), we note that confirmation of (138]) is consistent with our earlier 
discussion (see after (112jl ) that £ T is independent of coupling constant; in this dimensional 
window, temperature is a boundary effect in (imaginary) time and is the only temporal scale 
in the problem. 

In order to calculate a d , we rewrite the gap equation at criticality as 

-- 2 ^rjW|M^-l|, (39) 



T 2 ~ T 2 d J Q (2ny I u q " 2q 

V(d/2) 



where T d q d l dq iT d = fttt^t) is the d-dimensional volume measure. Rescaling A = a d T and 



q = uT, we obtain 

F d [a] = T 3 ~ d a 2 d h c (40) 

where 

, cothdVw 2 + oi 2 ) 1 



F d [a] = , =3 , , , - / u^dui \ 2 ~- ■ (41) 



For d < 3, the temperature prefactor on the right-hand side of ( [40]) vanishes T — > 0, so a 
consistent solution requires ad to satisfy 

F d [a d ] = 0. (42) 
At a small finite temperature, we can expand around a = ad + 5a (T), to obtain 

A(r)=ajr+ yy r " <43 » 
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Figure 5: Dependence of a<i(T — > 0) on dimensionality d. Inset: T-dependence of A/T for d = 2 
and couplings in the range 0.01 < 7c ( rf ~ 2 ) < 5.0; here Tn is the temperature scale where £ ~ a and 
we note that lim^^o a 2 is independent of j c . 

Thus in dimensions d < 3, the dominant low temperature behavior is independent of j c , 
the strength of the mode-mode coupling, which enters into the subleading temperature 
dependence. 

The necessity of separating out the singular part of equation (41) was pointed out to us 
by Chamati and Tonchev;— (41) was incorrectly treated in an earlier version of this paper. 
Following their approach, we can rewrite (41) as 

c " d ~ 1 r( d )r( 1 ~ d ) 



F d [a] 



{2^YT{d/2) 



oc 



ud _ ld ,coth(|V^)-l )+/ u ^ du 



Uo v v^ 2 + « 2 y Jo Vv^ + o 2 

(44) 

yielding 

fiW = — ^i— — r M ,- v roothtlvTO)-!, 12 ad _, ( r( ¥ < 



(2y/t) d T(d/2) J I v^T^ 1 J (2v^) d+1 ^(d-l)(3-d) 

(45)' 

provided (1 < d < 3). The first term in this expression is a smooth positive function of d 
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and a, whereas the second is a singular negative function of d with poles at d — 1 and d = 3. 
The numerical solution of -FrffaJ = can then be determined and is presented in Fig. [51 
We note that this result indicates that vanishes in the vicinity of d ~ 3 as aa ~ a/3 — d, 
consistent with previous e calculations.— 

In Figure 5 we display the dependence of ad on dimensionality 1 < d < 3. The 
temperature-dependence of the gap in two dimensions is shown in the inset of Fig. 5, 
where we see that lim<r^ a 2 — 0-96 is the same for all couplings. According to (]40j) and 
( )4T|) . we write a\ = 7 c F 3 [a] and solve for a 3 in the limit of upper cutoff u max = q m ax/T = 
2irT /T » « 3 , 

«3(r,7«)~y=g=, (46) 

where again we do not consider logarithmic corrections to a%. In the limit of strong coupling, 
«3 ~ [ln(-^^)l 1//2 is 7 C independent. For weak coupling, the situation relevant here, a s is 
indeed a function of 7 C but remains independent of temperature so that A ~ T according 
to (138j) : temperature-dependences derived here should therefore be in agreement with those 
found from a scaling perspective whenever direct comparison is possible. 



C. Temperature-Dependent Dielectric Susceptibility 

To provide an explicit illustration of the above calculations, we now use ( 1331) . and (135j) 
to numerically determine the temperature-dependent paraelectric susceptibility in the ap- 
proach to the quantum critical point (QCP) in d = 3. We obtain x _1 (T) = A 2 ~ T 2 for the 
approach r = r c in agreement with previous results and discussion. We note that a similar 
analysis in the vicinity of the classical phase transition leads to the expected Curie suscepti- 
bility (% _1 (T — > T c + >> 0) ~ T) since in this (high) temperature regime the Bose function 
in (136]) scales as — . We also remark that if we assume that to = ujq with no q-dependence 
then we recover the Barrett^ expression x^ 1 ~ A + 5coth^^; because the disperson is 
constant and q-independent this approach is not applicable near quantum criticality where 
the gap vanishes and the q-dependence becomes important. 

One more point needs to be considered before we proceed with our self-consistent Hartree 
theory. In the self-consistent Hartree theory (SCHT) of the ferro-electric phase, the polariza- 
tion field Pq acquires a nonzero value. Pq enters the Lagrangian Ce in ( 1221) as P = Pq + SP, 
where 5P are fluctuations of the polarization field around its mean value, P (P = in the 
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Figure 6: Temperature-dependence of the a) spectral gap and b) the dielectric susceptibility for 
three temperature scans defined in the schematic inset; here g = r — r c . 

paraelectric phase). The self energy (T2~4"l) becomes 

£ = 3 7c (P 2 > = 3 7c (P 2 + (<^ 2 » (47) 

as indicated diagrammatically in Figure 4. The equilibrium value Pq is easily obtained by 
introducing an electric field into the Lagrangian by replacing Ce — ► Ce + E ■ P, then seeking 
the stationary point 5S/5Pq = which gives (tPq + 37 c <5P 2 Po + 7c -Po 3 ~ E) = 0, or 

r + S - 2 7c P 2 = |- = (48) 

at zero electric field. According to (127|) . A 2 (r, T) = r + S(r, T), so that the spectral gap in 
the ferroelectric phase is 

Aj(r,T)=2 7c P 2 (r,T)>0. (49) 

In Fig. [6] a) we plot the calculated temperature-dependent spectral gap A(r, T) for three 
different values of r as indicated in its schematic inset. As expected, for (2) the spectral gap 
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Figure 7: T-g Phase Diagram as determined by a self-consistent analysis of the dielectric suscepti- 
bilty. The power law exponents are depicted in different colors via the function ^(T^—T y To ■ This 
expression selects the exponent 2 (blue region) for x~ l ~ T 2 (T c = for QCP), exponent 1 (green 
region) for classical Curie behavior x^ 1 ~ (T ~ ^c) an d exponent (red region) for a constant 
susceptibility. 

closes exactly at T = leading to a linear dispersion relation, u = q at the QCP. We note 
that in the quantum paraelectric (QPE), A (or x _1 ) is constant. In the quantum ferroelectric 
(QFE) again A is constant; though there exists a classical paraelectric-ferroelectric transition 
at T = T c where x" 1 ~ (T — T c ). The static dielectric susceptibility in the vicinity of the 
QCP (low T) is presented in the same three r regimes in Figj6] b) where we see that in 
the QPE regime x(T — > 0) saturates, at the QCP x(T) ~ T~ 2 and diverges as T — > 0. 
In the QFE, the susceptibility also saturates at low temperatures, though the Curie law is 
recovered in the vicinity of the classical transition at T = T c . 

Figure [7] shows the phase diagram that results from the self-consistent Hartree theory. 
This figure serves to emphasize how the strictly zero temperature QCP gives rise to a 
quadratic power law dependence of the inverse susceptibility on temperature over a sub- 
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stantial region of the T — g phase diagram. 

The crossover temperature, To, between Curie (x~ l ~ T) and Quantum Critical {x~ l ~ 
T 2 ) behavior in the susceptibility is defined by the expression 

r -£* < 5 °> 

where c^o = is the characteristic soft mode frequency, c s is the soft-mode velocity and a is 
the lattice spacing. Here we have assumed a simple bandstructure u(q) = ojq sin qa such that 
c s = ^j- = (uj a) cos gal so so that co> = — as stated above. The factor of 2ir in the 
denominator of (1501) results from the observation that the separation of the poles of the Bose 
and Fermi functions in the complex frequency plane is Ai/ n = 27r/csT, which sets the natural 
conversion factor between temperature and frequency to be 2irkB- Tq also corresponds to 
the temperature when the correlation length is comparable to the lattice constant (£ ~ a); 
here the correlation length £ = ~ g" 1 / 2 (see (j!7p ). Neutron scattering measurements^ 
of the dispersion relation indicate that the soft mode velocity in STO is c s ~ 10 4 m/s and 
the lattice constant has been measured^ to be asro — 3.9 x 10~ 10 m; therefore T ~ 30K. 
We note that with O 18 substitution, the ambient pressure Curie temperature^ is T c ~ 25K. 
Using the values of c s and asro from above, we get c^o ~ 2.6 x 10 13 Hz in SrTiO^. The typical 
frequency f2 = g 1 ^ 2 (spectral gap at zero temperature) at which one observes the change of 
behavior in the dielectric susceptibility (blue region) is thus from Figure 7, f2 ~ 10 -2 u^o = 
2.6 x 10 u Hz. Indeed, Raman scattering on ferroelectric SrTi 18 Os (T c = 25K) shows that 
the zero temperature Raman shifts is about 10 cm" 1 which translates into Qq ~ 3 x 10 n .£/z, 
in good agreement with our calculation. 



V. COUPLING TO LONG-WAVELENGTH ACOUSTIC MODES 



A. Overview 



In a conventional solid, broken translational symmetry leads to three acoustic Goldstone 
modes. At a ferroelectric QCP, these three modes are supplemented by one or more optical 
zero modes. This coexistence of acoustic and optic zero modes is a unique property of the 
ferroelectric QCP, and in this section we examine how their interaction influences observable 
properties. 
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The gap of the optical modes in a ferroelectric is sensitive to the dimensions of the unit 
cell and couples linearly to the strain field. This leads to an inevitable coupling between the 
critical optical mode and the long-wavelength acoustic phonons that must be considered. 
To address this issue, we consider the effect of a coupling 77 between the soft polarization 
and the strain field created by a single long-wavelength acoustic phonon mode. Softening 
of the polar transverse optic (TO) mode near the QCP enhances the effect of this coupling. 
Using dimensional analysis we find that the coupling between the TO and LA mode is 
marginally relevant in the physically important dimension d = 3, and thus can not be 
ignored. The main result of the analysis is that the acoustic phonons act to soften and reduce 
the quartic interaction between the optic phonons. Beyond a certain threshold 77 > rj c , this 
interaction becomes attractive, leading to the development of a reentrant paraelectric phase 
at finite temperatures. We note that such a coupling to acoustic phonons has been considered 
previously— and here we are rederiving and extending prior results in a contemporary 
framework. 

B. Lagrangian and Dimensional Analysis 

We introduce the coupling of the polarization (P(x, r)) and the acoustic phonon r)) 
fields as a coupling of the polarization to strain — 7/V0P 2 ; we then write the Lagrangian^ 

as 



where £e[P] is our previous Lagrangian without acoustic coupling given in (1221 . Here 
the constant 77 is the coupling strength to the acoustic phonon; the latter's dynamics are 
introduced in the bracketed terms of (l5Tl) . Since we are using units in which the velocity of 
the soft optical phonon is one, c = ^ is the ratio of the acoustic to the soft optical phonon 
velocities. We will discuss the restoration of dimensional constants in (151j) when we make 
comparison to experiment in Section V. F. 

We begin with a dimensional analysis of the couplings to assess their relative impor- 
tance in the physically important dimension d = 3. In order to do so, we introduce the 
renormalization group (RG) flow by rescaling length, time, momentum and frequency 



C E [P, 0] = C E [P] + g [{9r<f>f + 5 2 (V0) 2 ] - 77V0P 2 



(51) 




T = 



j, q' = <?A, 



v 1 = z/A, 



(52) 
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with constant A > 1 representing flow away from the infrared (IR) limit of the QCP, that 
is flow from small to large momentum and frequency. In terms of the rescaled variables x' 
and t', the action ([3]) with Lagrangian (1511) in d + 1 dimensions becomes 

S[P,(f>] = J drj d d xC E [P,4>] 

/^/ A r r1 r i 

= / dr' d d x'A d+1 {^A- 2 [(d T ,P) 2 + (VP) 2 + {d T 4f + (V'c0) 2 

+ \^lP 2 + \lcP A -r/A-V0P 2 }. (53) 

We emphasize that we write = r — r c as the coefficient of the P 2 term in the Lagrangian 
Le[P] (|22|) . entering (15T1) in (1531) . since our RG flow starts from the QCP (r = r c ). Rescaling 
P, <p, Qq, 7 c and n, so that the action (|53l) assumes its initial form, we write 

P' = PAV 0' = M^, (^)' = ^A 2 , 7 ^ = 7cA 3 - d , v' = V^~^, (54) 
which leads to 

/•/3/A /• 

S[P,(t>] = J dr' J d d x'C E [P',(j)']. (55) 

Now the fields, the mass term and the coupling constants flow to new values leaving the 
action unperturbed. We remark that the upper cuttoff in the imaginary time dimension is 
replaced by infinity as the temperature T ~ ^ approaches zero. 

Analyzing the RG expressions in (154|) . we find that the term grows as we flow away 
from the QCP IR limit; therefore it is a relevant perturbation parameter independent of 
dimension d. This is consistent with the fact that Qq = r — r c = g tunes the system away 
from the QCP. Similarly we find that couplings 7 C and rj grow (relevant) in dimension d < 3, 
decrease (irrelevant) in dimension d > 3, and don't change (marginally relevant) in d = 3. 
We see that in this case (d = 3) the coupling to acoustic phonons (7/) is equally important 
as the mode-mode coupling (7^) and thus has to be included to the Gaussian model. 

Let us now briefly summarize what we know about 7 C before we proceed to the discussion 
of the acoustic coupling rj. In section IV B we found that the spectral gap A is independent 
of 7 C for dimensions 1 < d < 3 in the zero temperature limit (see Fig. 5). This is in 
agreement with the above results, where 7 C is a relevant perturbative parameter; more 
precise analysis^ shows 7 C flowing to the nontrivial Wilson-Fisher fixed point 7*. Here all 
the system properties become functions of 7* + S r y c ~ 7*, and so are 7 c -independent. On the 
other hand, in dimensions d > 3 and d — 3, j c flows to zero (with logarithmic corrections in 
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Figure 8: Diagrammatic representation of the self-energy that includes coupling to both optical 
and acoustic phonons. Here j c is the renormalized coupling, including the exchange of an acoustic 
phonon. 

the marginal case). In these cases the system properties are functions of S^c and thus are 
7 c -dependent; we have already seen an example of this behavior in the specific case of the 
d = 3 spectral gap. 

C. Gap Equation 

We are now ready to explore how the system's low-temperature behavior changes in the 
presence of acoustic phonons in dimension d = 3. Let us look first at the LA phonon field 
<p. Following the procedure of Section IV A, we find the acoustic Green's function and 
dispersion relation from (151]) to be 



We emphasize the P 2 -dependency of the new interaction term, — 7/V0P 2 , in the La- 
grangian (15TT) . Therefore it contributes to the polarization self-energy as an additional term 
inside the brackets of (1231) . This new contribution arises due to nonzero second-order per- 
turbation and is schematically shown in Figure 8, where the solid line represents the soft 
polarization TO Green's function (1251) and the dashed line represents the LA Green's func- 
tion (1561) . We note that the interaction represented by a dot in the self-energy consists of a 
contribution each from the coupling 7 C and t]. Thus we can write the polarization self-energy 
S as a sum of these two terms 



D(q) = D(q,iv n ) = [{iv n ) 2 



~2 21 — 

c q ] 



(56) 



v a (q) = cq. 



(57) 



E(r,T) = E 7c (r,T) + E,(r,r) 
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n ^ 

+ VT^ J ^ q "G(q,iu n )D(q,iu n ), (58) 

where £ 7c is the Hartree self-energy (1301 previously calculated in Section IV A. We remark 
that the q 2 term in the integral for H v arises due to form of the interaction (V0). Converting 
the Matsubara summation to a contour integral, deformed around the poles z p = ±u p (q) 
and z a = ±u a (q) in the dispersion relations of the polarization (I26I) and acoustic phonon 
( 157|) modes respectively, we can rewrite in the forrn^ 

E, (r ,T) = -4, y |____ + ___). (59 ) 

At the quantum critical point, where r = r c and T = 0, the dispersion cj p (g) = q and 
^s(^p) — n B(u a ) = so that 

^°>- 4 ' 2 /(S^T15' (60) 

Using (|28|) . we write the gap function (as in IV A) as 

A 2 = nl + Al + Al, 

d 3 q ({n B {u p (q)) + \] 1 



" 7c " ,D J (2vr) 3 I u p 2q 

A 2 = _ 4 „2 f _&q_ n 2 ( + , [n B Mq)) + l] _ 

' " M (2tt)» 9 ^ i^K-afl ^K-^] 25[5+l]g3 I 

where A 2 c has been already defined in (1531) . 

We emphasize that the 7 C and rj terms in (l6Tj) have opposite signs in their contribu- 
tion to the spectral gap A. The negative coefficient of rf reflects the fact that it emerges 
from second-order perturbation theory; physically it is due to thermally enhanced virtual 
excitations caused by coupling between polarization TO and LA phonon modes. 



D. Deep in the Quantum Paraelectric Phase 

Let us first explore the effect of the acoustic coupling 77 deep in the QPE region of the 
phase diagram (see inset of Figure 9). Here g » and A >> T « 0. In this regime, we 
write 

x ~ 1 = A 2 = nl + D(A)-A( v )^ (62) 
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Figure 9: Schematic temperature-dependence of the static dielectric susceptibility where coupling 
to a long-wavelength acoustic phonon is included in the calculation; inset indicates phase trajectory 
and region of corrections due to acoustic coupling deep in the QPE phase (yellow). 

with 

A(rj) = -y- J j^un B (cu), (63) 

where derivations of A(rj) and -D(A) are presented in Appendix A; for our purposes, the 
key point is to note that limA^o-D(A) = 0. Setting A(rj) ~ rf = 0, we recover a con- 
stant expression for \ as a function of temperature in the QPE phase, consistent with our 
previous derivations from Section IV. For rj ^ 0, the dielectric susceptibility acquires new 
low-temperature behavior. The quartic temperature term in (1621) . —A(rj) ^ , drives the 
inverse susceptibility at low temperatures; such a "bump" in the susceptibility (or "well" in 
the inverse susceptibility, see Fig. 9) due to acoustic phonon coupling has been considered 
previously^. It is then natural to enquire whether a finite r\ could eventually drive the 
inverse susceptibility to zero (or negative) values. Here we show that this is not the case. 
We start by looking for a solution of ( 1621) with x" 1 = A 2 = 0, and show that such a solution 
does not exist. Indeed at rf = 0, x _1 i n the QPE phase is nonzero as we saw in Section 
IV. At rf 7^ 0, growth of last term in (|62j) exceeds all bounds and cannot equate a constant 
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term fi 2 , (notice that D(A)| A =o = °)- The 

inverse susceptibility therefore remains positive 

deep in the QPE phase with Xmin ^ 0- 

We note that when the temperature increases so that A ~ T and we are no longer 
in the QPE phase (red in Fig. 7), we enter the "tornado" region of the QCP influence 
(blue in Fig. 7) where x~ l ~ A 2 ~ T 2 , as was shown in Section IV. At this point the 
quadratic temperature-dependence dominates and coupling to the acoustic phonons becomes 
negligible; as a result a turn-over in the inverse susceptibility from — T 4 to +T 2 -dependence 
occurs (see Fig. 9). 

E. Quantum Critical Temperature-Dependent Dielectric Susceptibility 

We already know that there exists a classical phase transition at T c for g < and rj = 0; 
for 7] ^ could this line of transitions enter the g > part of the phase diagram and result 
in a reentrant quantum ferroelectric phase near the g = QCP? In order to explore this 
possibility, we study the temperature-dependent susceptibility near the QCP (at g = 0) 
and find that unstable behavior is possible. Next we follow the line of transitions, where 
X^ 1 = A 2 = and show that its behavior is changed for rj > r/ c . 

We begin with x{T) in the vicinity in the quantum critical regime where g = (trajectory 
2 in Figure 6); here = g = and g~co>~T>0at low temperatures. Taking 77 = 0, the 
spectral gap ( ISTl) becomes 

K = ^j dqqn B { q /T) = a lc T* = (64) 

and we recover the quadratic temperature dependence, x^l = ^ 7c ~ T 2 , that was derived 
in Section IV B. 

With 7] ^ 0, the 77 contribution to the gap becomes 

^-^S d «^ (»»(*^ - ^f 3 ) - <«) 

For both cases c ^ 1, the expression under the integral in fl65l) is positive (see Appendix B 
for specifics), which results in a negative coefficient for A 2 . Adding both 7 C and 77 terms in 
the gap equation flBTl) . we write the expression for the dielectric susceptibility 

X' 1 = A 2 = (a 7c - ^ 2 )T 2 = (h- f3 V 2 )T 2 (66) 
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Figure 10: The transition line T c {g) for different values of rj, the acoustic coupling constant; for 
r) > r/ c a reentrant quantum ferroelectric (FE) phase emerges. The phase boundaries result from 
numerical solution of the gap equation (A(rj ^ 0) = 0) discussed in the text; the parameters used 
here are: 7 C = l,c = 0.9, r\ c = 0.62 and {g > i] c ,i] < g c } = {0.63,0.59}. 

where a and (3 are explicitly calculated in Appendix C. When the coefficient of T 2 is zero, 
namely when 




the phase boundary line (x 1 = 0) becomes vertical in the approach to the QCP; when 
X] > r] c , it "meanders" to the right leading to reentrant behavior. 

F. Details of the Phase Boundary ( x _1 = 0) 

We now follow the phase transition line, denned by x" 1 = (A = 0) out to finite 
temperatures. From Section IV we know that there is a classical ferroelectric-paraelectric 
phase transtion at g < at Curie temperature T c (g); it is depicted as a solid line in Fig. 6, 
where the dielectric susceptibility diverges, x = A -2 — > oo. Our results in Section IV are 
for rj = 0, and we study the effect of rj > r] c > on this transition line. 

To do this, we look for a solution to the gap equation ( 16T|) . when A(T C , rj) = 0, which 
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yields the transition line T c (g,r]). When the spectral gap closes, the dispersion relations of 
the TO soft polarization and the LA acoustic modes both become linear (oo p (q) = q and 
u) a (q) = cq). Inserting these values into (IBTj) and setting A = 0, we obtain 

Arf f qmax r . n B (cq/T c 



2vr 2 ^ = 37 c J dqqn B (q/T c ) - _ J dqq{n B (q/T c ) - 



c 



{/•Umax PUmax , Tig (cu) ~) 1 
37 c y duun B (u) - — — — J duu{n B (u) z — jj (68) 

for the equation determining T c (g). At low temperatures, we note that we recover the scaling 
relation Q% = g ~ T 2 since both integrals become proportional to T 2 {u max = >> 1). 
At high temperatures n B (u) ~ -, so the r.h.s. of (16"Sj) becomes proportional to T c , and we 
recover the classical behavior g ~ T c . 

Figure 10 shows the T c (g) transition line. For rj > rj c w 0.6, the transition line "wanders" 
into the g > region, leading to a reentrant quantum ferroelectric phase. Such reentrance 
suggests the possibility of nearby coexistence and a line of first-order transitions ending in 
a tricritical point, but the confirmation of this phase behavior requires a calculation beyond 
that presented here and will be the topic of future work. 

In order to make direct comparison with experiment, we must now restore dimensions 
to our coupling constant and more generally to our Lagrangian (15T|) . We start by explicitly 
restoring all physical coefficients to the Lagrangian, as follows 

(3F = [ d ^L 



2 



h 

{d T Pf + c 2 (VP) 2 ] + r -fP> + ^P 4 - MV0)P 2 + P- 



(69) 



where c s and c a are the soft optical and acoustic phonon velocities respectively and where 
P and <j) are the un-rescaled physical polarization and phonon displacement fields. Then by 
writing 

-=x, c ^p^ = p\ ^4> 2 = ^ 2 (70) 

Co h h 



we obtain floTj) . the rescaled Lagrangian, 
(3F = J d 3 idrC E (P,(l)) 

c e (p, <p) = \ [{d T pf + (VP) 2 ] + r -p 2 + |p 4 - ?KV0)P 2 + \ [{m 2 + £ 2 (V0) 2 ] (71) 
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where 



c 2 s r D h l h _ c a 

jL = C E , r = -, 7 = ^7D, ^7^VV D ' c =- ( 72 ) 



In the dimensionless units used in this Section, we found that 




(73) 

where a = ~ and (3 = y J^f^ ) • Using (1721) . we can now rewrite this critical coupling in 
dimensionful terms as follows 

[P "5/2 

V n 



h 3 \ 4 Vc 2 + c+l 




For SrTiOs, the acoustio^ 2 - 1 ^ and the soft-mode^ velocities have been measured to be c a 
8000m/s and c s ~ lOOOOm/s respectively so that c = 0.8; the crystal mass density is 
5.13g/cm 3 = 5.13 x 10 3 kg/m 3 . The value of 7 C has been measurecU^ 1 ^ using ferroelectric 
Arrott plots of E/P vs P 2 to be eo7 C D = 0.2 m 4 /C 2 . Inputting all these numbers and 
rj c = 0.6 into our dimensionful expression for r) C D, we obtain 

tjSto = 5 74 x 10 io Jm/C 2 (75) 

as the dimensionful critical coupling to be compared with experiment. 

Next we estimate the experimental value of rj in SrTiO^ as^ 4 - r) STO ~ ® where Q and s 
are the typical magnitudes of the electrostrictive constants and the elastic compliances 4 ^ 4 ^ 
respectively; here we use the values^ 4 - Q = 0.05 ^ and s = 3x 10~ 12 Thereofore we 

obtain 

VsTO = 1.7 x 10 10 Jm/C 2 (76) 

so that from our analysis we observe that rjsro < VcD° f° r the SrTiO^ system. However, 
there are two points of uncertainty here that we should emphasize: (i) we use experimental 
values for SrTi 1& 0% as they are not yet available for 5 , rTi 18 3 ; (ii) we use values of Q and s at 
room temperature, and these quantities need to be determined at low temperatures. Despite 

28 



the roughness of our estimate, it is reasonable to assume r\ is not changed dramatically by the 
issues raised in (i) and (ii). We encourage further experimental investigations of SrTi 18 Oz 
at low temperatures to clarify this situation. 

G. Translational-Invariance as Protection against Damping Effects and Singular 
Interactions 

Our analysis of the effects of acoustic coupling has been limited to a Hartree treatment 
of the leading self energy. This approach neglects two physical effects: 

• Damping, the process by which a soft mode phonon can decay by the emission of an 
acoustic phonon 

• The possibility of singular interactions induced by the exchange of acoustic phonons 

Similar issues are of great importance in magnetic quantum phase transitions in metals, 
where the coupling of the magnetization to the particle-hole continuum of electrons intro- 
duces damping.— For example, in the simplest Hertz-Moriya treatment of a ferromag- 
netic quantum critical point, damping by the electron gas gives rise to a quadratic Lagrangian 



where the term linear in \u\ is a consequence of damping by the particle- hole continuum. 
This term plays a vital role in the quantum critical behavior; by comparing the dimensions 
of the q 2 term with the damping term, we see that [u] = q 3 , which means that the temporal 
dimension scales as z = 3 spatial dimensions under the renormalization group. This has the 
effect of pushing the upper critical dimension down from 4 — 1 = 3 to 4 — 2 = 1 dimensions. 
In addition to this effect, the coupling to the electron- hole continuum also introduces non- 
local interactions between the magnetization modes, casting doubt on the mapping to a 4 
field theory. 

Fortunately, translational invariance protects the ferroelectric against these difficulties. 
Translational invariance guarantees that the soft mode can not couple directly to the dis- 
placement of the lattice; instead it couples to the strain, the gradient of the displacement, 
according to the interaction Hj = —r]V<f)P 2 . When we integrate out the acoustic phonons, 



of the form 




q,v 



(77) 
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the induced interaction between the soft-mode phonons takes the form 



,2 



V(q, v) = -V I - , (78) 



where the numerator result from the coupling to the strain, rather than the displacement. 
The presence of the q 2 term in the numerator removes the "Coulomb- like" 1/q 2 divergence 
at small q, protecting the soft mode interactions from the development of a singular long 
range component. 

A similar effect takes place with the damping. To see this, we need to examine the 
imaginary part of the self-energy appearing in the Gaussian contribution to the action, (1231) . 



q,v 



v 2 + q 2 + r + E v (q) 



|P(g»| 2 . (79) 



Damping results from the imaginary part of self energy, S"(g, uS). To compute the damping, 
we generalize T, v given in (1591) to finite frequency, obtaining 

E v (q,z) = 2 V j 7 —^k\ u _ ^ — + — J7 - ^ — + ( z ++- z )^ (80) 



(2ir) 3 lu p [(z - Up) 2 - uj 2 ] u a [(z - u a ) 2 - u. 



where we have used the short-hand u a = u a (k), uj v = u p (q — k), n a = UB(uJ a (k)), n p = 
fiB{^p{q — k)). The imaginary part of this expression at zero temperature, for positive z/, is 
then given by 

/d 3 k k 2 
TTT^ 5(u-U) a -u) p ). (81) 
{2-K) A UJ p UJ a 

We can determine the small q, uj behavior of this damping rate by simple dimensional 
analysis. The dimension of the right-hand side is [g 5 ]/[co> 2 ] ~ q 2 , so the damping rate must 
have the form 

ImS^^-r/ 2 ^^,^, (82) 

where a more careful analysis of the integral reveals that F j^j is not singular at either 
small momentum or frequency. The most important aspect of this result is that the scattering 
phase space grows quadratically with frequency and momentum, so that it does not dominate 
over the other terms in the action (1791) . The scaling dimension of frequency remains the 
same as that of momentum, and thus the upper-critical spatial dimension remains as d — 3. 
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VI. DISCUSSION 



A. Logs, Dipolar Interactions and the Barrett Formula 

Before summarizing our results, let us briefly touch on a number of topics closely related 
to our work which we have not yet discused; more specifically they include logarithmic 
corrections in the upper critical dimension, dipolar interactions and the use of the Barrett 
formula for quantum paraelectrics. As we have already noted in Section VB, the polarization 
mode-mode interaction j c , and coupling to the acoustic phonons 77, are both marginally 
relevant in the dimension of physical interest d = 3. Thus logarithmic corrections to the 
scaling relations (III) have to be included; we have already seen their appearance in the 
expression for in (146]) . The correction to scaling of the free energy near the classical 
ferro-paraelectric phase transition in four dimensions is^ 

/dft 7c) = /oft 7c) [1 + ^cHU/t)] 1 / 3 , (83) 

where t = \ T ^ Tc | is the reduced temperature, /oft 7c) = ( \t/t^/ 2 ) * s ^ e sc& l m S form 
of the free energy with a universal scaling function $, to is the reduced Debye temperature 
for the soft mode (1501) and 7c is the polarization mode-mode coupling at QCP. Since x — 
f^|fi=o, we have 

X = Xo[l + 9lcln(t /t)f 3 , (84) 

where Xo ~ t -1 . By applying the quantum-classical analogy (III), we write at the upper 
critical dimension, = 3 (d + z — 4; z — 1), 

/ g m(<?,7c) = ZoG?,7 c )[l + 9 7 cMV4)] 1/3 , (85) 

where g = uj^ is the Debye frequency for the soft mode squared, fo(g, 7 C ) has the same form 
as before, and g is the tuning parameter. By setting x = q^i\e=q, the dielectric susceptibility 
becomes 

X = Xo[l + 9 7c ^(Vs)] 1/3 , (86) 

where Xo ~ 9 1 ~ T~ 2 . The temperature-dependence of x with logarithmic corrections is 
then found by making the subsitution g ~ T 2 in (186j) . and these results are identical to those 
found previously using diagrammatic techniques^. An analogous procedure can be used to 
find the logarithmic corrections to other thermodynamic quantitites. 
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We note that here we assume the upper critical (spatial) dimension d" = 3; however 
if we include uniaxial dipole-dipole interactions, we will have d u c = 2. Basically this is 
because when all dipoles point in the (same) z-direction, the TO polarization frequency (126]) 
becomes^ 

W 2 (g) = g 2 + A 2 + /4, (87) 

where (3 is a constant, and we derive (I87j) in Appendix D. We note that the last term of 
(187)) is specific to the uniaxial (e.g. tetragonal) case and is not present for isotropic dipolar 
interactions. Applying simple scaling, we obtain 

where the constants b, b k > 1 represent flow to the infrared (IR) limit of the QCP. We show 
in Appendix D that in order for ( 187)) and ( 1881) to be satisfied simultaneously, k must equal 
2 so that q z "counts" for effectively two dimensions [d s ^ e — d + 1), so that for a quantum 
uniaxial ferroelectric the total effective dimension is d e ff = d s ^ e + z = (d+1) + z = d + 2 
with d" = 2 since then we obtain d e ff = 4. 

At this time, it is not known whether SrTi l% 0% is cubic or tetragonal at low temperatures. 
In any case, we expect the samples under study to be structurally multi-domain so that 
averaging over long length-scales will make them effectively cubic; thus uniaxial dipolar 
interactions do not need to be considered. The observed T 2 behavior of x m the vicinity of 
the QCP supports this contention (i.e. d s e p ^ ce = 3); for d s ^ a p = 4, a different T-dependence 
{x~ x ~ T 3 ) is expected 2 ^ for a QPE so that a reexamination of the underlying model would 
be necessary to match experiment. Until details of the samples are known, this situation 
cannot be ascertained. We note that such T 2 dependence of the inverse susceptibility has also 
been observed 2 ^ in mixed crystal ferroelectrics KTa\- x Nb x O^ and Kai- y Na y TaO^ where 
uniaxial dipolar interactions are not important, and we encourage further low-temperature 
studies of these systems. 

A consistent discrepancy between the observed low-temperature dielectric susceptibility 
and the Barrett formula^ has been observed in the quantum paraelectric phase.— >^ Here 
we emphasize that the discrepancy occurs when the system gets very close to the QCP; thus 
it provides a measure of the tuning distance to the QCP. Because the optical polarization 
mode softens as the system approaches the QCP, with the gap vanishing completely here, 
the momentum dependence in the dispersion relation (126)) becomes important. It is exactly 
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for this reason that the Barrett formula, that assumes a constant dispersion relation, u = uq, 
breaks down close to the QCP. 

The Barrett formula^ works well deep in the QPE phase (V D), where the gap is much 
bigger than temperature. One such example is KTaO% (KTO), which remains paraelectric 
down to zero temperature, but in contrast to SrTiO^ (STO) shows a much lower value of the 
zero temperature dielectric susceptibility (xkto ~ 4000, Xsto ~ 24000)^^. The closer the 
system is tuned to the QCP, the smaller is the spectral gap and the bigger is the dielectric 
susceptibility Therefore, STO sits much closer to the QCP than KTO, and indeed KTO 
shows a nice fit to the Barrett formula^. Notice that by plugging d> into ( 1331) . we get the 
Barrett expression, 



where Xi = uq, and M and To are fitting constants. 
B. Summary and Open Questions 

Let us now summarize the main results of the paper. Here our aim has been to character- 
ize the finite-temperature properties of a material close to its quantum ferroelectric critical 
point; we have rederived and extended previous theoretical results using scaling methods 
and self-consistent Hartree theory. In the process we have made an analogy between tem- 
perature as a boundary effect in time and the Casimir effect, and have used this to shed 
light on both problems. Using simple finite-size scaling, we have presented straightforward 
derivations of finite-temperature observables for direct comparison with experiment, and our 
approach has yielded a scaling form xi 00 ) = ^-^XfO which serves as an additional probe of 
To, the soft-mode Debye temperature-scale where we expect crossover between Curie (T) 
and Quantum Critical (T 2 ) behavior in We emphasize that this scaling method is useful 
in this system where z is low (z = 1); otherwise if z is higher, the system is usually well 
above its upper critical dimension where this approach is inappropriate. Next we've used 
self-consistent Hartree methods to determine the T — g phase diagram and the crossover 
between classical and quantum behavior. In particular we see the influence of the quantum 
critical point on the susceptibility at finite temperatures, and we can put in materials pa- 




(89) 
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rameters to determine the size of its basin of attraction. Finally we include coupling to an 
acoustic phonon and find that it affects the transition line; for such couplings greater than 
a threshhold strength there is a reentrant quantum ferroelectric phase. 

Naturally these results suggest a number of open questions and here we list a few: 

• The presence of a reentrant phase suggests the possibility of nearby phase coexistence, 
a tricritical point and a line of first order transitions. This is a particularly appealing 
scenario given that recent experiments^- suggest coexistence of QPE and QFE in 
SrTi 18 Os and is a topic we plan to pursue shortly. 

• If indeed there is a tricritical point and a line of first-order phase transitions, could 
there also be a metaelectric critical point in the g — E plane analogous to the metam- 
agnetic situational^ in some metallic systems? There is indication that an analogous 
metaelectric critical point occurs in a multiferroic system,— so this is a question driven 
by recent experiment. 

• What happens when we add spins to a system near its quantum ferroelectric critical 
point? Would the resulting multiferroic have particularly distinctive properties? 

• Similarly what type of behavior do we expect if we dope a quantum parelectric in the 
vicinity of a QCP? There is by now an extensive body evidence that electronically 
mediated superconductivity is driven by the vicinity to a magnetic quantum critical 
point, phenomenon of "avoided criticality" , whereby superconductivity in the vicinity 
of a naked magnetic quantum critical point£^£. In such systems, the metallic transport 
properties develop strange metallic properties that have been termed "non-Fermi liquid 
behavior"-^. This raises the important question, as to what, if any, is the ferroelectric 
counterpart to this behavior? In particular - how does the presence of a soft mode 
affect the semi-metallic properties of a doped quantum critical ferro-electric, and does 
a doped ferroelectric quantum critical point also develop superconductivity via the 
mechanism of avoided criticality? 

We believe that we have only begun to explore the rich physics associated with the 
quantum ferroelectric critical point, a simple setting for studying many issues associated 
with quantum criticality that emerge in much more complex materials. Furthermore the 
possibility of detailed interplay between theory and experiment is very encouraging. 
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VIII. APPENDIX A: D(A) AND A(rj) 

We derive expressions for D(A) and A (77) (J63~j) using the gap equation flSTj) deep in the 
QPE region (D), where g » and A >> T « 0. Collecting all "-"-terms under integrals 
of A 2 c and A 2 in ( IBTi) . we obtain the expression for D(A), 

1 1 



D(A) 




ui a [u; 2 - u 2 } q 3 c[l - c 2 } 
_ 37c T V 2 , V 2 T 

— T~2 1 ~ ^~2 i2 2 i3 ' 



h = 
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'0 



A 2 + q 2 [l-d 2 ] q 2 [l-c 2 } 



m a :c 



/3=/ d " l ' [ WT7v^¥^\ ~ ^i^T 1 ' (90> 







Notice that \im.^oD(A) = 0, since all three integrals Ii, I2 and ^3 become zero at zero gap. 
We split the integrals Ii (i = 1,2,3) into two parts, Ii = J Q nA + J^^T" 1 , where nA » A. Since 
q» A in the second integral part, we neglect its A dependence and get a zero contribution. 
Thus, only the first integral part contributes, and D(A) becomes a function of A only, with 
no temperature dependence. 

Next we show that the second Bose-Einstein ns(o; a (g)) term under the integral of A 2 in 
(|6T|) results in the form A{rj) in equation (163]) . 
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f d 3 u „ T 4 

J Wf U ~ ~ A{V) A2 ' (91) 



where u = q/T. Notice that we approximate A 2 >> q 2 (l — c 2 ) in the second line of (l9Tj) . For 
low momenta, this is indeed the case. For large momenta, g>>A>>T^0, we neglect A 
in (I9~T|) and the integral becomes 

A V 2 



2tt 2 5(1 - 5 2 ' 



J dqqn B (cq). (92) 



In the limit q >> T, n B (cq) ~ e~ cq ^ T and fl92|) becomes exponentially small (~ T 2 e~ cq ^ T ) 
and can be neglected. Similarly, we neglect the rest of the terms in the gap function (l6Tj) 
with Bose-Einstein thermal distribution n B (u p (q)). Deep in the QPE phase A >> T, so 
that n B {uj p {q)) ~ e~ A / T at low momenta, or n B (to p (q)) « e -iiarg e /T & ^ j ar g e momenta. In 
both cases A,q [arge >> T, the integrals containing n B (ujp(q)) become exponentially small 
and so are negligible. 

IX. APPENDIX B: INTEGRAL (JBS5 IS POSITIVE FOR c ^ 1 

We also show that the expression under the integral in (I65I) is positive for the two cases, 
c 0s 1. First, assuming that c < 1, cq < q (positive g's) and n B {cq/T) > n B (q/T) we write 

{ nB(5/r) _!!^} i > (i-I )ns (c,/r) » 

I c J c z — 1 c <r — 1 

n B {cq/T) > 0, (93) 



£(£+1) 

which we note is positive. Similarly, for c > 1, cq > q and n B {cq/T) < n B (q/T), we write 

; M , /r) _ } i > (1 _ i) ns(9/r) i 

L C J c z — 1 c & — 1 

n B (g/T) > 0. (94) 



c(c+ 1 

which is also positive. Therefore the integral in (161)1) is positive in both cases. 

X. APPENDIX C: a AND ARE CONSTANTS 

To evaluate the quantities a and /3 in (I64p and (I65p . we make a change of variables to 
u = q/T, and u = cq/T respectively. The expressions for these two constants then become 

a = —r / duun B {u) = - 

2n Jo 4 
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where we have taken the limits of integration to infinity and used the result 
lo 
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XL APPENDIX D: DIPOLE-DIPOLE INTERACTIONS IN UNIAXIAL FERRO- 
ELECTRICS 

The interaction energy between two dipoles pt and pj siting on two sites t\ and r} respec- 
tively is 

3 47re |f| 3 ' W 

where n is a unit vector in the direction of the vector r = rj- — r*j. From ( l96i) . we find the 
total dipole-dipole interaction potential to be 



5 n h 3r a r 6 



(97) 



i j, a, 6 

where r = |r|, and a, b label vector coordinates. After we perform a Fourier transform, the 
interaction potential becomes 



^ = rE^- (98) 



a 



where g = \q\ refers to the momentum-dependence of W(q). Assuming that all dipoles point 
in the sanies-direction in the uniaxial case, we find that the dipole potential 

W(q) ~ & (99) 

W(q) contributes to Lagrangian (l5Tj) . Le[P, 3>] — > Le[P, 3>] + W, so that the TO polarization 
frequency fl26|) then reads^ 



^ 2 (g) = c 2 g 2 + A 2 + /4, (100) 

g z 



where we introduce constant of proportionality (3. 

We show that (IH71) and (EEL conditon k = 2. Let us assume that > 1. Then 
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Since we also rescale frequency uo{q) (j57|) by a constant, expressions, q 2 and are to be 
proportional. This leads then to the condition 

k = 2. (102) 
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